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Abstract 

A spatially confined stability analysis is reported for the cylinder wake at Reynolds numbers 
190 and 260. The two three-dimensional instabilities at these Reynolds numbers are shown to be 
driven by the flow just behind the cylinder. 
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Consider the von Karman vortex street generated by flow past an infinitely long circular 
cylinder. This flow is timeperiodic and two dimensional for Reynolds numbers between 
approximately 47 and 189 . (The Reynolds number is Re = U^d/v, where L/^ is 

the fluid velocity far from the cylinder, d is the cylinder diameter, and v is the kinematic 
viscosity.) At Reynolds number 189 the 2D vortex street becomes three dimensionally 
unstable S fl fl fl Q. A nntnerical stability analysis of the flew np to ReynoHs n Um be r 
300 has determined two separate bands of linearly unstable modes The first appears 
at Re = 189 with a spanwise wavelength of 4 cylinder diameters and the second appears 
at Re = 260 with a spanwise wavelength of 0.8 diameters. These linear instabilities are 
related to 3D shedding modes known as mode A and mode B, first observed experimentally 
by Williamson 

In this paper I revisit the stability analysis of the vortex street. I show that small 
regions of the full flow just behind the cylinder are responsible for the 3D linear instabilities 
despite the fact that the actual linear modes extend many cylinder diameters downstream 
of the cylinder. This is important for two separate reasons. The first is that this limits 
the regions which should be analyzed either to clarify instability mechanisms or to suppress 
the instabilities if so desired. The second reason is that the stability analysis of small flow 
regions is computationally very efficient compared with a stability analysis of the full flow 
field. Thus this approach provides a method for quickly obtaining approximate stability 
information about a complex flow. 

The computational approach is described fully elsewhere [^J . I summarize briefly the main 
points with focus on the important aspects for the current work. First direct simulations 
of the incompressible Navier-Stokes equations are used to obtain 2D, time-periodic wake 
flows: XJ(x, y, t + T) = XJ(x, y, t) where T is the wake period. A spectral-element method 
is employed on a computational domain shown in Fig. ^ The length scale is the cylinder 
diameter. The size of the computational domain depends on the Reynolds number and the 
accuracy required The domain shown gives sufficiently accurate results for the purposes 
of this work. 

The boundary conditions imposed for the computation of U are as follows. On the top, 
bottom, and left boundaries uniform flow is imposed U = (U, V) = (1,0). At the right edge 
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FIG. 1: Spectral-element computational mesh. The base flow is computed on the full domain 
and stability computations are performed on the full domain and on a variety of subdomains. The 
subdomain outlined in bold and shown as an enlargement has inflow length Li = 1.5, outflow 
length L a = 4, and cross-stream height L y = 1.5. 

a (Neumann) outflow boundary condition for the velocity and pressure P is used: 



No-slip conditions are imposed at the cylinder surface. 

The next step is a Floquet stability analysis of 3D disturbances to the 2D wake. This 
analysis is based on the evolution of infinitesimal 3D perturbations u'(x, y, z, t) of the 2D flow 
XJ(x, y, t). The equation for such a perturbation is obtained by linearizing the Navier-Stokes 
equations about U: 



where p' is the perturbation to the pressure that enforces V ■ u' = 0. 

In this work I shall primarily consider Eq. (J2J) posed on subdomains such as that illustrated 
in Fig. ^ These subdomains are characterized by their inflow and outflow lengths Lj and 
L , measured with respect to the cylinder center, and their cross-flow range — L y < y < L y . 
The base flow U appearing in Eq. (J2J) is obtained by restricting the base flow from the full 
domain to the appropriate subdomain. 



dU/dx = 0, dV/dx = 0, P = 0. 



(1) 





Boundary conditions on u' are required by Eq. (J2J. Consider the standard case in which 
u' is computed on the same domain as U. Then the requirement that U + u' satisfy the same 
boundary conditions as U gives homogeneous boundary conditions on u': u' = (0, 0, 0) on 
all domain boundaries except the outflow boundary where u' satisfies the analogy of Eq. ([[J. 
One can view these boundary conditions as the requirement that the perturbation be zero 
at the inflow and lateral boundaries of the domain but that perturbations can advect out 
of the domain. Based on this reasoning we impose these homogeneous boundary conditions 
when computing u' on subdomains. We are interested in examining perturbations u' which 
are zero on all boundaries of the subdomain except the outflow boundary. Perturbations are 
allowed to advect out of the subdomain. 

From the linearized equations the spectrum of Floquet multipliers fi and corresponding 
Floquet modes u can be found ^j]. Exponentially growing perturbations correspond to 
multipliers outside the unit circle in the complex plane > 1). Because the geometry is 
homogeneous in the spanwise direction, Floquet modes decouple in spanwise wavenumber (3. 
The Floquet multipliers are thus computed as a function of (3 on a number of subdomains for 
Re = 190 and Re = 260, just above the critical Reynolds numbers for the 3D instabilities. 

Figure |2] summarizes the study in terms of multiplier spectra while Figs. 00 shows rep- 
resentative Floquet modes. A large number of subdomains have been studied but only 
representative cases are shown near the minimum dimensions necessary to capture the 3D 
instabilities. 

Figure El primarily demonstrates the effect of outflow length L Q on the Floquet multiplier 
spectrum. The inflow length and cross-stream height are fixed at minimal values of Li = 
and L y = 1.5. One sees that L Q = 3 is sufficient to capture the destabilizing Floquet mode. 
Increasing L a further has minimal effect. For Re = 190, decreasing L D to 2.25 causes a 
substantial change to the spectrum - the subdomain is evidently too small to capture all 
the relevant flow features which drive the instability. In the case of Re = 260, decreasing 
L D to 2.25 causes only a small change to the spectrum and the subdomain is essentially still 
large enough to capture relevant flow features. Decreasing the outflow length to L = 1.5 
causes the spectrum in both cases to deviate wildly from the correct behavior (at Re = 260 
the multipliers become complex with large magnitude and cannot be plotted on the scale of 
Fig. |2J). Figs. 00 shows just how well the Floquet modes are captured by the subdomain 
computations even though the Floquet mode on the full domain extends far downstream of 
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FIG. 2: Summary of Floquet multipliers at Re = 190 and Re = 260 (as labeled) for j3 near the 
critical value in each case. Solid curves are from computations on the full domain. Points are from 
subdomains with L, = 0, L y = 1.5, and L a = 3 (squares), L Q = 2.25 (crosses), and L a = 1.5 (dots). 
The multipliers for Re = 260, L Q = 1.5 are off the scale of the figure. 

the cylinder P|. 

Increasing the inflow length Li has some effect on the spectrum of multipliers in the 
case of Re = 190. However, the effect is very small at the critical wavenumber where the 
multiplier branch has largest magnitude. Increasing the inflow length has almost no effect 
in the case Re = 260. Increasing the cross-stream height has very little effect in either case. 
It is possible, in fact, that the cross-stream height could be reduced slightly from that used 
here without compromising the accuracy of the multipliers. 

It is reasonable to conclude that the physical region responsible for the instability is 
localized within the subdomain in those cases where the Floquet multipliers and Floquet 
modes from subdomains are close to those from the full domain. This localization of the 
instability is similar to the transition from steady flow to vortex shedding which takes place 
at the primary wake instability. In that case there is a local, absolutely unstable region 
extending approximately three to four diameters downstream j^, very similar to the 
scales observed here. In fact the current computational analysis is similar in spirit to that of 
Hannemann and Oertel :8] and can be thought of as situated between the local analysis of a 
parallel wake profile and the global analysis previously reported for the secondary instability 
to three dimensionality j^. 

Finally, it is worth considering the stability of the downstream portion of the flow. In 
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FIG. 3: Floquet modes at Re = 190 with (3 = 1.6. The main figure shows (a portion of) the mode 
computed on the full domain. Insets show modes computed on the subdomains indicated. (The 
inflow boundary Lj = and the lateral boundaries L y = 1.5 are the same for all subdomains.) All 
modes are plotted as spanwise vorticity contours at the same fixed time. The spanwise vorticity 
contours ±0.4 of the base flow are also plotted. 

this case the base flow is projected onto the subdomain with 2.25 < x < 25 and —4 < y < 4. 
Both at Re = 190 and Re = 260 the flow in this downstream region is found to be very 
stable. The maximum modulus of the Floquet multipliers are ~ 0.28 at Re = 190, 
(3 = 1.6 and ~ 0.19 at Re = 260, (3 = 7.5. 

In this paper I have shown how a local Floquet stability analysis of small regions of the 
wake just behind a cylinder are sufficient to capture both 3D instabilities observed in wake 



transition. The specific results are of interest 



physical mechanisms of 3D transition 



or future attempts to understand better the 



0, H n, 12, Q 



because the results set limits to 



the regions which could drive the instabilities. The general method is of interest for future 
computational studies of flow instabilities. 

Much of this work was performed while visiting the Courant Institute of Mathematical 
Sciences, NYU. I thank members of the Institute for their hospitality. 
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FIG. 4: Floquet modes at Re = 260 with /? = 7.7 using the same conventions as Fig. El 
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